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Abstract 

As an alternative to the paradigm of slow roll inflation, we propose an 
extended scenario of the matter bounce cosmology in which the Universe 
has experienced a quasi-matter contracting phase with a variable background 
equation of state parameter. This extended matter bounce scenario can be 
realized by considering a single scalar field evolving along an approximately 
exponential potential. Our result reveals that the rolling of the scalar field in 
general leads to a running behavior on the spectral index of primordial cos¬ 
mological perturbations and a negative running can be realized in this model. 
We constrain the corresponding parameter space by using the newly released 
Planck data. To apply this scenario, we revisit bouncing cosmologies within 
the context of modified gravity theories, in particular, the holonomy correc¬ 
ted loop quantum cosmology and teleparallel F(T) gravity. A gravitational 
process of reheating is presented in such a matter bounce scenario to demon¬ 
strate the condition of satisfying current observations. We also comment on 
several unresolved issues that often appear in matter bounce models. 

Pacs numbers: 04.50.Kd, 98.80.Bp, 98.80.Jk 


1 Introduction 

The Matter Bounce Scenario (MBS) (see E El II SlU for recent reviews on a 
variety of bouncing cosmologies) suggests that the Big Bang singularity is replaced 
by a nonsingular bounce. It is essentially characterized by the Universe being 
nearly matter dominated at very early times in the contracting phase (to obtain 
an approximately scale invariant power spectrum) and evolving towards a bounce 
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where all the parts of the Universe become in causal contact (6), solving the horizon 
problem, to enter afterwards into a regular expanding phase, in which it matches 
the behavior of the standard hot big bang cosmology. As an alternative to the slow 
roll inflationary paradigm, the MBS can be free from the problem of the initial 
singularity Q that inflation models suffer from. Also, the potential for the scalar 
field in MBS does not need to be extremely flat as what is required in inflationary 
cosmology HI. 

In order to obtain a viable MBS model that can compete with the inflationary 
paradigm, it is expected that the underlying model can satisfy a variety of theoret¬ 
ical and observational constraints. Moreover, there exist several conceptual issues 
that are not clear in the frame of the MBS. We list these points in the following and 
discuss the conditions of model building in the MBS in the present work. 

First, today’s cosmological measurements, such as the Planck data released in 
2013 (Planck2013) S |T0l as well as in 2015 (Planck2015) ifTTl H2ll . have pre¬ 
cisely determined the amplitude of the power spectrum for primordial curvature 
perturbations to be V g = 2.2 x 10 -9 . This amplitude, in bounce models, is often 
associated with the energy scale of the bounce as well as the process of primordial 
perturbations evolving through the nonsingular bouncing phase llT3l . Accordingly, 
the theoretical result of the power spectrum of primordial curvature perturbation 
calculated from any bounce model has to match with the observational data. 

Second, according to the Planck2013 data, at ler confidence level (C.L.) the 
spectral index for curvature perturbation and its running, namely, n s and a s , are 
constrained to be 0.9603T0.0073 and — 0.0134T0.009, respectively ©[10]]. These 
results are obtained upon the standard Lambda Cold Dark Matter (ACDM) model 
with a parameterized power spectrum and thus, are general to any inflation and 
alternative models. It is well-known that the ways to obtain a nearly scale invariant 
power spectrum of primordial perturbations are either a quasi de Sitter phase in the 
expanding phase or a nearly matter dominated phase in the contracting phase fl4l . 
However, in an exact MBS with the background equation of state being w = 0, 
one finds n s = 1 for the single field model and therefore, conflicts with observa¬ 
tions (which is similar to the fact that an exact de Sitter inflation was ruled out by 
the data). In order to improve the scenario to correctly explain observations, as 
in inflationary cosmology where a dynamically slow roll period was introduced, 
one expects a quasi matter dominated phase in bouncing cosmology with the back¬ 
ground equation of state deviating from zero by the condition w = ^ -Cl, where 
P and p the pressure and the energy density of the Universe, must be considered at 
early times in the contracting phase. 

Thud, it is important to take into account the observational constraint of prim¬ 
ordial non-gaussianity upon early Universe models. So far there is no evidence 
pointing to the existence of these nonlinear fluctuations |[T5l fl6ll . As a result a 
large number of inflation models are ruled out by this observational fact. Thus, the 
no-detection of primordial non-gaussianity is expected to tightly constrain bounce 
models. For instance, it was studied in detail in lfl7ft that primordial non-gaussianity 
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generated in the simplest version of MBS is not sizable enough to be probed in 
cosmological observations. However, for the bounce models achieved by includ¬ 
ing a spatial curvature term, the amplitude of the nonlinear fluctuations could be 
dangerous to explain latest observations llTSll . In general, the signature of the prim¬ 
ordial non-gaussianity via the calculation of the 3-point function remains unclear 
in a generic picture of MBS since the calculation relies on the detailed process for 
primordial fluctuations evolving through the bouncing phase which contains uncer¬ 
tainties depending on different mechanisms of generating a nonsingular bounce. 

Fourth, the latest CMB experiments including the BICEP2/Keck Array and 
Planck data have constrained the tensor-to-scalar ratio to be r < 0.12 with a pivot 
scale of 0.05 MpcU 1 at 2 a C.L. |fT9| . When applied to inflation models, due to a 

slow roll consistency relation r = 16e ll20l (where e = —jyi — \ \ VJ ' s ^ 1C 
slow roll parameter), it requires the inflaton’s potential to be very flat and hence 
indicates a fine tuning issue EO. In the simplest model of MBS, the amplitude of 
tensor fluctuation is comparable to that of curvature perturbation and thus the value 
of r is too large ll22ll . Accordingly, one may expect certain dynamical mechanisms 
to be implemented in the MBS in order to enhance the amplitude of curvature 
perturbation to be consistent with the data, namely, a curvaton mechanism ff23l . 
Also, it is possible to consider the effects of modified gravity theories such as from 
loop quantum cosmology to depress gravitational waves during the bounce If24l . 
Therefore, the no-detection of primordial gravitational waves can also impose a 
bound on various bounce models. 

Fifth, a Universe whose background dynamics is realized by a primordial scalar 
field (not by isotropic fluids) has to reheat via decaying into light particles that will 
thermalize to match with the standard hot big bang expansion. Reheating could 
be produced due to the gravitational particle creation in an expanding Universe 
ll25l [26l . In this case, an abrupt phase transition (a non-adiabatic transition) is 
needed in order to obtain sufficient particle creation that thermalizes, producing 
a reheating temperature that fits well with current observations. In bouncing cos¬ 
mologies, the gravitational particle creation is natural to be implemented since the 
Universe would have experienced several phases including contracting, bouncing 
and expanding ones ll27l . and then, it is necessary to examine whether the reheating 
temperature is compatible with the current data Il28ll . 

Sixth, bouncing cosmologies often suffer from a dangerous growth of primor¬ 
dial anisotropy of which the effective energy density scales as a -6 in contracting 
phase. This is known as the famous Belinsky-Khalatnikov-Lifshitz (BKL) instabil¬ 
ity ||29l . A solution to this problem can be realized in the ekpyrotic scenario f30l . 
in which the energy density of the dominant matter scales with a~ q with q 6. 
To keep the scale invariance of primordial curvature perturbation, a model of MBS 
including a period of ekpyrotic contraction was recently suggested in OTTl and the 
evolution of primordial anisotropy was analyzed in lf32l . However, it is important 
to be aware of this issue in other mechanisms of nonsingular bounces, such as tak¬ 
ing into account nonlinear matter contribution to smooth out the anisotropies li33il 
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or by including higher order curvature term If34l . 

Seventh, studies of distant type la supernovae lf35ll reveal that our Universe 
is expanding in an accelerating way. A viable cosmological model needs to ac¬ 
commodate with this late time acceleration, which usually is incorporated with a 
cosmological constant, or by a scalar field (see If36l for recent reviews). There are 
other ways to implement the current cosmic acceleration, for example using F(R) 
gravity (see for instance Il37ll38l0 . Then it becomes important to question whether 
a bouncing solution can be realized in these different models ll39ll . 

From the above considerations, the aim of the present work is to address, from 
a critical viewpoint, part of those points. In particular', dealing with a single scalar 
field, we provide a clear definition of the quasi-matter domination regime, which 
is the key to obtaining a running behavior of the spectral index, and we will see 
that this period is characterized essentially by two quasi-matter contraction para¬ 
meters. We perform a detailed calculation of the power spectrum of primordial 
perturbations, which allows us to obtain the expression of the spectral index and 
its running in terms of these two parameters associated with quasi-matter contrac¬ 
tion. Once we have obtained these expressions, we can investigate the parameter 
spaces of some models, showing that these models of MBS could fit well with latest 
cosmological observations. We also present a detailed study of the gravitational re¬ 
heating in the MBS when the background is depicted by holonomy corrected Loop 
Quantum Cosmology (LQC). More precisely, we show that when one considers a 
massless field nearly conformally coupled with gravity, the reheating temperature 
provided by this model is compatible with current observations. Afterwards, we 
comment on several unclear issues that remain in specific models of MBS, such as 
nonlinearities produced during the bouncing phase, the Jeans instability in the case 
of holonomy corrected LQC, as well the non-local Lorentz invariance of the scalar 
torsion in bouncing teleparallel models. 

The units used throughout the paper are: h = c = 8 irG = 1. 


2 The phase of quasi-matter contraction 


In order to obtain a nearly matter dominated contracting phase, a simple way is to 
consider a potential of the form 

V{<p) = U 0 e~^ 3(1+ “ )M , (1) 


where cr is a very small parameter. Then, the exact solution to the Friedmann and 
conservation equations, 


¥>(*) = 


1 

y / 3(l + a) 


In 


r3(l + a) 2 
-2(1-a) 



( 2 ) 


depicts a quasi-matter dominated Universe with equation of state P = ap. As we 
will see in next section, this background corresponds to a spectral index equal to 
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n s — 1 = 12a with vanishing running. 

To depict a viable phase of quasi-matter contraction, we consider, for a flat 
Friedmann-Lemaitre-Robertson-Walker (FLRW) geometry, the Friedmann and con¬ 
servation equations for the homogeneous paid of a single scalar field, 

H ' 2 = 1(y + V )’ <P + 3H<p + Vv = 0. (3) 


Assuming that there was a quasi-matter domination in the contracting phase, i.e., 
ip 2 = 2V ==> <p = V v , the above background equations become 


f H 2 = \V f U 2 = la 2 V 

\3H<p + 2V ip = 0 ^ \ 3Hp' + 2a% = 0, 


( 4 ) 


where T~L is the conformal Hubble parameter and a prime denotes a derivative with 
respect to conformal time. 

This stage can be depicted by the equation of state parameter. 


P 

w = — 
P 


2 /1 

3\2 + WJ 



V) 


-1, 


( 5 ) 


where P and p are the pressure and the energy density, which arc related to the 
spectral index, and characterize this regime through the condition that w <C 1. 

Since a potential of the form Voe~^^ generates a phase of exact matter- 
domination, we would like to reexpress our potential V, for negative values of the 
field, as 

V(<p) = Voe^a+PM) } (6) 


and thus we obtain 


w “ 2 ((pf{ip)) v , 


(V) 


where we have assumed that |/(<^)| <C 1 and \ipfip(<p)\ <C 1. Further, it is conveni¬ 
ent to introduce another parameter as follows, 


S 2 = ^ = ~ = ~V3{vf{^))<pcp =. ( 8 ) 

which is found to be related to the running of the spectral index in next section. 

To demonstrate the small deviation from an exact matter-domination in the 
above model, we consider examples by choosing the functions f(ip) = - 

(/3 > 0) and f(tp) = ^ ri , for negatives values of the fields satisfying \tp\ 2> 1. 
Then we obtain 


w 

\/3 ip ’ 


w = -2(2n — 1) 




( 9 ) 
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and 



( 10 ) 


respectively. 

Moreover, in analogue with inflationary cosmology, we can also define an ef¬ 
fective number of e-folds before the end of the quasi-matter contracting phase as: 
a(N ) = e N af, where a/ is the value of the scale factor at the end of this stage. 
With this definition, in the quasi-matter approximation, the number of e-folds is 
calculated as 



( 11 ) 


which, in terms of the function f(ip), can be expressed as, 



( 12 ) 


As a final remark to this section, we note that it is well-known that when one 
only considers a single canonical scalar field in spatially flat cosmologies, the frame 
of General Relativity (GR) forbids bounces from the contracting to the expanding 
phase, as is best seen by looking at the Raychaudhury equation in the flat FLRW 
space-time H = — \<p 2 < 0: the Hubble parameter always decreases, so it is 
absolutely impossible to pass from negative to positive values. For this reason, 
when the matter part of the Lagrangian is given in terms of a single scalar field, 
one is led to consider a non-canonical scalar field involving Horndeski operators 
such as model constructions in lf3Tll32l . or to use cosmologies beyond the realm of 
GR, namely, holonomy corrected LQC, teleparallel F(T) gravity, or Fill) gravity. 

3 Current MBS with a single scalar field: Background 
Review 

In this section we review the main MBS with a single scalar field in the frames 
of LQC and teleparallel F(T) gravity, where the modified Friedmann equation is 
given by an ellipse in the plane (H, p) that depicts a Universe moving clockwise 
along it |6!|- 

3.1 Holonomy corrected LQC 

The main idea of LQC is that it assumes a discrete nature of space which leads to 
consider, at quantum level, a Hilbert space where quantum states are represented by 
almost periodic functions of the dynamical part of the connection l40l . However, 
the connection variable does not correspond to a well defined quantum operator 
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in such a Hilbert space and therefore the re-expression of the gravitational part of 
the Hamiltonian in terms of almost periodic functions (holonomies) is needed. It 
might be executed from a process of regularization. This new regularized Hamilto¬ 
nian introduces a quadratic modification (p 2 ) in the Friedmann equation at high 
energies HU . which can give rise to a nonsingular bounce when the energy density 
becomes equal to a critical value bellow the Planck energy density. More precisely, 
the holonomy corrected Friedmann equation in LQC, which depicts, in the plane 
(H, p ), the ellipse given by 


h 2= I{^£)’ <13) 

where p c is the so-called “critical density” (the energy at which the Universe 
bounces). 

3.2 Teleparallelism 

It was also noticed in iflZl that a nonsingular bouncing solution can be derived in 
theories of teleparallel gravity. Teleparallel theories are based on the Weitzenbock 
space-time. This space is M 4 , with a Lorentz metric, in which a global, orthonor¬ 
mal basis of its tangent bundle given by four vector fields {e*} has been selected, 
that is, they satisfy g(ei,ej ) = rpj with // = diag (—1,1,1,1). The Weitzenbock 
connection V is defined by imposing that the basis vectors e, be absolutely parallel, 
i.e., Ve; = 0. 

The Weitzenbock connection is compatible with the metric g, and it has zero 
curvature because of the global parallel transport defined by the basis {e*}. The 
information of the Weitzenbock connection is carried by its torsion, and its basic 
invariant is the scalar torsion T. The connection, and its torsion, depend on the 
choice of orthonormal basis {e^}, but if one adopts the flat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric and selects as orthonormal basis {eo = do . e\ = 
-<9i,e2 = -c? 2 , e 3 = -c? 3 j, then the scalar torsion is T = —6 H 2 , where H = - 
is the Hubble parameter, and this identity is invariant with respect to local Lorentz 
transformations that only depend on time, i.e. of the form e) = A ^{t)ek (see 

MM). 

With the above choice of orthonormal fields, the Lagrangian of the F(T) theory 
of gravity is expressed as: 


C t = V{F(T) + C m ) , (14) 

where, once again, V = a 3 is the element of volume, and Cm is the matter Lag¬ 
rangian density. Accordingly, the Hamiltonian of the system is given by 

Ut = ~ F(T) + p) V , (15) 
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where p is the energy density. Imposing the Hamiltonian constrain Ft = 0 leads 
to the modified Friedmann equation 

p = -2^^r + F(T) = G(T) , (16) 

which defines a curve in the plane ( H , p) by applying the relation T = —6 H 2 . It 
is interesting to note that Eq. ( |T6j ) may be inverted, and therefore, each curve of 
the form p = G{T) defines an F(T) theory through the following expression: 


F(T) 


\f—T 

2 


G(T) 

Tsf^T 


dT. 


(17) 


In order to produce a cyclically evolving Universe, let us take the F(T) theory 
arising from the ellipse that defines the holonomy corrected Friedmann equation in 
LQC (Eq. (p~3|>). To obtain a parametrization of the form p = G(T), the curve has 
to be split in two branches 

e = G±{T) = 08) 

where the branch p = G-jT) corresponds to 77 < 0 and p = G + (T) is the branch 
with H > 0. Applying Eq. ( | 1 7[ ) to these branches yields (lf6l) 

F±{T) = ± \j — r ~^~ arcsin + G ±( T ) ■ ( 19 ) 


3.3 Phase space dynamics 

In holonomy corrected LQC and teleparallel F(T) gravity with a single scalar field, 
namely ip, the dynamical equation of the background is given by the conservation 
equation 


ip + 3F[±ip + V v (<p) — 0 , 


( 20 ) 


where H = - x /f (1 - £ 


= t/f(l-£)i"the 


;J , ^ p , in the contracting phase, // — y , 

expanding one. 

The goal of this Section is to prove that, when one deals with a potential that 
defines a quasi-matter domination at early times in the contrac ting phase, i.e., with 
the form V(<p) = VoeF^vCL+fM)^ a ll the solutions of Eq. (20) depict, at early 
times, a matter dominated Universe. 


To show this property, first of all, note that Eq. (20 1 has a solution with the 
following asymptotic behavior at early times 


p{t) 



when t —> —oo, 


( 21 ) 













because one can check that — In y §Voij is solution of the equation 

</j + 3iT±<^> + V v (ip) = 0 , (22) 

with H- = —\J\ and V(tp) = V 0 e^. 

Secondly, at early and late times one can disregard holonomy corrections, and 
we can reduce the analysis to large values of \ip\, and thus, take the potential 
V(tp) = Then, performing the change of variable ;p = — \nip, the 

corresponding non-perturbed Klein-Gordon equation (or equivalently, the conser¬ 
vation equation) reads 


with 


dip 

dtp 


F±W, 


(23) 


F±W = '^(U 2 + V 0 

4-0 \d 




(24) 


where in F ±, the sign + (resp. —) means that the Universe is in the expanding (resp. 
contracting) phase. Note that, Eq. ( [23] ) describes two (one for the contracting and 
other one for the expanding phase) one-dimensional first order autonomous dy¬ 
namical systems, that are completely understood calculating critical points. These 

critical points are ip + = 


| Vo f° r the expanding phase and ip- = — y f Vo for the 

contracting one, meaning that in the contracting phase ip- is a global repeller, i.e., 
at very early time all the solutions have the asymptotic behavior ip(t) —> ip- when 

t —> —oo. 


Finally, the point ip- corresponds to the solution ip = — ^ In ( — y |Vq t 


which satisfies p = ip 2, = ^, what proves that all the solution of Eq. (|20| depict, 


at early times, a matter-dominated the Universe in the contracting phase. 


4 Detailed calculation of the power spectrum in quasi¬ 
matter bounce 

Common to all these cases is the Mukhanov-Sasaki B31 equation for scalar per¬ 
turbations, in Fourier space, which can be expressed as 

v 'k + 

where for very low energy densities and curvatures, c 2 = 1 and z = a fj = a^. 
The explicit expressions for z and c 2 in the cases of F(T) and Fiji) gravities were 
derived in li46ll and li47l . respectively. In contrast, for the holonomy corrected LQC, 


Ik 2 - Z - ) v k = 0 , 


(25) 
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the form of 2 is the same as in GR, i.e., z = a^-. However, the expression for the 
square of the sound speed parameter differs from unity at high energies |i48l . As 
we will discuss later, in the holonomy corrected LQC, the value of the square of 
the sound speed parameter becomes negative during the super-inflationary phase 
H > 0, which happens when the energy density is greater than one half the critical 
density. 

During the quasi-matter dominated contraction, which occurs in the contracting 
phase at low energies and curvatures, one has = — (a detailed derivation of this 
relation was done in Il49l0 . which means 


— ^ H' + U 2 

z 


(26) 


Accordingly, solving Eq. ^ for a constant w, i.e., taking H = |(1 — 3 w) and 
inserting into d26b, we get 


z" 2 

- = -2(l-9«;), 

Z Tj z 


(27) 


under the assumption of \w\ <C 1. We refer to lf50l for a more generic analysis. 

It is clear from this result that the Mukhanov-Sasaki equation (25 1 during the 
quasi-matter dominated contraction can be approximately expressed as 


// 1 

v k + 


k 2 - f v 2 - - 

r/ z 


Vk = 0 ; v=--6w. 


(28) 


Assuming the initial conditions of primordial perturbations to be vacuum fluctu¬ 
ations, one then obtains the solution to Eq. 


Vk = 




For modes well outside of the Hubble radius A:| 77 I <C 1, (28) becomes 

l/o 1' 


Vk ~^y - 4) Vk = 0 ’ 

and then the dominant mode of the solution is given by 

v k = c 2 (k) \n\*~ v . 


(29) 


(30) 


(31) 


On the other hand, to calculate the power spectrum we will choose a pivot 
scale, namely A;*, and let 77 * be the time at which the pivot scale crosses the Hubble 
radius in the contracting phase. Then we could write the scale factor as follows 

HA V 



a(rf) = a. 
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2 


(32) 










where a* and H * are the values of the scale factor and Hubble parameter, respect¬ 
ively, at the crossing time. The approximation on the r.h.s. comes from the fact 
that in the quasi-matter dominated contracting phase, we have aH = |. Hence, 

since in the quasi-matter domination one has z(rj) = \/3 JY^w)a(r]) = y/3a(r]), 
the solution ( [37| ) can be written as follows, 

v k = 3\Z3-^-\rj*\~^~ u C 2 (k) (z(r)) j ■ (33) 


For modes well outside of the Hubble radius, the solution ( |29l ) should match 
Eq. ( |33| ). Using the small argument approximation in the Hankel function and 
expression we find 


£k(v) 



1 

16 



f v dr) 

\H*\ J-oo Z 2 (fj) 



(34) 


These modes will re-enter the Hubble radius at late times in the expanding phase, 
when the Universe will be matter dominated. Then, the power spectrum is given 
by 


V^k) 


27 kl f n dfj \ 2 fk\ 3 ~ 2u 
64vr 2 i7 2 Vd-oo - 2 (.V)) 


( 35 ) 


Provided that the Universe is always dominated by the same matter field through¬ 
out the contracting and expanding phases, and assuming that the usual procedures 
apply all through the non-singular bouncing phase, which is justified due to the 
analycity of £k(v) ( see Eq. (34 1 ) in the bouncing phase, one can evaluate the above 
equation and eventually obtain the final formula for the power spectrum corres¬ 
ponding to scalar perturbations at the pivot scale A:*: 


V^K) 


27 fcg / f +°° d v \ 2 
647T 2 i7 2 \ 7 _oo z 2 (y)J 


(36) 


at late times of cosmic evolution. 

The calculation of the power spectrum, in general, can be performed numer¬ 
ically as pointed out in ll24l . However, dealing with LQC or teleparallel F(T) 
gravity, in the simple case of a matter dominated Universe throughout the whole 

background evolution, that is, for the scale factor a(t) = a* (^H*t 2 + , 

where p c is once again the critical density, i.e., the value of the energy density at 
the bouncing time, the calculation could be done analytically giving as a result 

Vt(K) = ^ — , (37) 

9 p p i 

in the case of holonomy corrected LQC iBTl . and 

'Pdk*) = ir—C 2 , (38) 

9 p p i 
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where C = 0.9159 is Catalan’s constant and p p i is the Planck density, for the case 
of teleparallel F(T) gravity | [46l . 

It is important to note that in the above calculation we have assumed that the 
Universe is dominated by a single matter field throughout the whole primordial 
epoch. However, a much more realistic Universe has to involve other matter com¬ 
ponents, such as radiation and a cosmological constant l39l . In addition, one often 
introduce an ekpyrotic matter field to dilute unnecessary anisotropies generated in 
contracting phase 11521 . In these cases, the amplitude of primordial curvature per¬ 
turbation will be different and it is possible to obtain the correct amplitude which 
satisfies the CMB observations even when p c is of the order of the Planck energy 
density ||39l . 

In order for the theoretical results provided by both models (Eqs. (37) and ( |38| )) 
to match with the CMB data = 2 x 10 -9 , the value of the energy density at the 
bouncing point has to be at least of the order p c ~ 10 ~ 9 p p i or even larger. 

From Eq. ([35]) and following the definition of the spectral index for scalar 


perturbations, one can derive the following expression: 


n, - i = 


^f) ^— 2 ,=i2„. 

dink 


(39) 


Note that one can further expand the spectral index by involving a running index 
via the parametrization, n s = n s (k*)+a s (k*) In (jr 'j ■ The running of the spectral 
index, a s , is then obtained as 


a, = 


dn q 


2 n’ 


24 w’ 


(1 - 3m) = —485 2 , (40) 


d\nk) k=aH (l + 3m)7f 
where we have used Eqs. ([5]) and ([8]). 

Remark 4.1. The quasi-matter domination was first addressed in H39I in the con¬ 
text of the ACDM model. In this context 


w = -—- < 0, (41) 

PA + PCDM 

where p,\ is the energy density corresponding to the cosmological constant A and 
PCDM I s die energy density of the Cold Dark Matter (CDM). When the CDM dom¬ 
inates the evolution of the Universe, the value of w is negative and close to zero 
and it is easy to show that the spectral index is given by n s — 1 = 12 w. In this 
specific case, after applying the approximate relations u/ = 3 TL{w + w 2 ) = 3 TLw 
and ppz = —\ (3 w + 2), one can calculate the running of the spectral index and 
finds 


OL s 


( dn s \ 
\dhtk) 


k=aH 


12 w"H 

W 


-72 w > 0 . 


( 42 ) 


As a result, the ACDM bounce model predicts a positive running behavior l\39\l . 
which is not ruled out yet by observations. Note also that the same happens if one 
mixes CDM with an isotropic fluid with negative pressure P = wp (w < 0). 
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In the same way, for tensor perturbations one has ^ = 'H'+'H : = ^ (1 — 9m), 
and thus, the following power spectrum is obtained: 


V T (k) 


32 kl ( r + °° dri \ 2 (k_\ 3 ~ 2u 

81lT 2 H 2 \J _oo Zt(v)J \ k *J 


(43) 


where, for very low energy densities and curvatures, zt = a. The exact expres¬ 
sion of zt in holonomy corrected LQC was obtained in ll53l . in teleparallel F(T ) 
gravity in Il46l . and in modified F(R) gravity in 1541 . 

In the MBS, the spectral index for tensor perturbation and its running are 


ut = 


dP T {k) 

din k 


= n s - 1, 


OIT — 


driT \ 

dhlk )k=aH 


(44) 


and the ratio of tensor-to-scalar perturbations, is given by 


r = 7M*0 8/ g^)M 2 

V d k ) 3 V/-^ z~ 2 (ri)dri) 


5 Fitting the spectral index and the running parameters 


The Planck2013 data released in 2013 constrain the values of the spectral index 
and its running to be 


n. 


= 0.9603 ± 0.0073 , = -0.0134 ± 0.0090 , 


at lcr |(9]. In models of slow roll inflation (such as monomial, natural, hilltop 
and plateau potentials) and in the MBS with a quasi-matter domination, in general, 
1 — n s is of order N~ 1 (N being the number of e-folds before the end of the corres¬ 
ponding regime), while the running parameter is of order N~ 2 and, consequently, 
one has a s ~ (1 — n s ) 2 . This situation is compatible with the Planck2013 data. For 
example, in the case of the MBS, for large negative values of ip, with the potential 
V(tp) = e^ v (—ip)P (i.e., for f(ip) = one h as 


n, - 1 = 


24 p 
\/3 ip' 


Ota - 


48 p 

ip 2 - 


(46) 


which means that a s = —4^(1 — n s ) 2 . This situation is compatible with a large 


range for /3 > 0. Choosing for instance B = j]- the constrains (461 are compatible 
with the Planck2013 data for 


p E (-17.81,-12.28). 

Dealing with the MBS model f(p) = , as has been introduced in the first 

section, one has 


n, 


24 

~ 1 = -2 ’ 

if 2 


OLa — 


96^3 


( p° 


(47) 
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which leads to the relation a 2 = 2{1 — n s ) 3 . A simple calculation shows that (47) 
fits well the Planck2013 data for values of the field in the range 


(p e (-27.24,-22.60) . 


(48) 


V 0 


The same happens, in MBS, with exponential potentials of the form V = 
~~—with B > 0. In this case one has 

l+e — p w\ 1 


n< 


_ 1 - 

Vs ’ 


a s = -48 B 2 e~^\ 


(49) 


which lead to the relation a s = 2\/3/3(n s — 1). This relation is compatible with 
the Planck2013 data for a wide range of values of /3, for example, when B = — 

On the other hand, in slow roll inflation one also has to take into account the 
latest CMB constrain r < 0.12, and the fact that more than 50 e-folds are needed 
to solve the horizon and flatness problems of standard big bang cosmology. By 
studying the (n s , r) plane, the Planck2013 data is able to constrain the parameter 
space of inflation models very well, preferring potentials with a concave shape 
(V w < o) Col. Moreover, these models could be further constrained if the data of 
the running parameter a s were included 1551 . 

Effectively, in slow roll inflation, (see for instance lf56l for a review of infla¬ 
tionary cosmology) the commonly used parameters are: 


e 


H 2 2 V V J ’ 


V u 


fj = 2e - — “ 4JW 
1 2 He V 


and the spectral index and its running are given by 

n s — 1 = 2f) — 6e, a s = 16ery — 24e 2 — 2f, 


(50) 


(51) 


where we have introduced the second order slow roll parameter 


£" = 



rj = 


4^? ip ip 

V 2 ' 


(52) 


Moreover, in inflationary cosmology, the tensor-to-scalar ratio r, is related to the 
slow roll parameter e, in the way 


r = 16e, (53) 

which does not happen in the MBS. 

Dealing with the running, a simple calculation leads to the relation 

a s = ^(n s - l)r + ^r 2 - 2£ . (54) 

Then, choosing for instance, the ACDM+r+a:. s model from Planck2013 combined 
with WP and BAO data, which gives the following results n s = 0.9607 ± 0.0063, 
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r < 0.25 at 95% C.L. and a s = -0.021±°;°] 2 (see table 5 of [10]). We will 
consider n s at 2a C.L. and take the conservative bound r < 0.32 (see figure 4 of 
mi Since the minimum of the function \ (n s — l)r + ^r 2 is reached at the point 
(n s = 0.9481, r = 0.1384) given a value greater than —0.0018, one obtains the 
bound 


> -0.0018 - 2£, (55) 

meaning that plateau potentials such as V(ip) = Vo ^1 — ^ (Hill-Top Inflation 

(HTP)) J57l, V(<p) = V 0 (l - (Double-Well Inflation (DWI)) El, V(<p) = 

Vo ^1 + cos (Natural Inflation (NI)) Ii59ll . or V(ip) = Vo ^1 + a In ^cos ) 
(Pseudo Natural Inflation (PSNI)) li60l . when one considers values of the running at 
la C.L., are disfavoured by Planck data because for all of them £ < 0. In fact, the 
distance from the theoretical value of the runnig to its expected observable value, 
namely 2 ), is larger than 1 . 6 cr. 

A distance larger than 1.6<r is also obtained for the potential that leads to Ex¬ 
ponential SUSY Inflation (ESI) V(ip) = Vq (1 — e~ pip ) 11611 . for Power Law In¬ 
flation (PLI) whose potential is given by V(ip) = V§e~ w | [62l . for Kaller Moduli 
Inflation I (KMII) K3l given by the potential V(ip) = Vo (1 — aipe~^), in Witten- 
O'Raifeartaigh Inflation (WRI) ll64l with potential V(<p) = Vo In 2 ^ j, and for 

general hill-top potentials such as: V((p ) = Vo ^1 — ^ (Small Field Inflation 

(SFI)) 1651 and V(tp) = Vo ^1 + (^j ^ (Valley Hybrid Inflation (VHI)) Il 66 l with 
p > 3, because in these cases one has £ < i] 2 , that is, 

p / 3 \ i 3 r 2 

> - l(n 8 - 1) + —rj - ~(n s - l ) 2 = > -0.0018, (56) 

where we have evaluated a s , as a function of n s and r at the absolute minimum, 
namely n s = 1 — ^ with r = 0.1384, in the rectangle [0.9481, 0.9733] x [0,0.32]. 

For potentials such as: V(<p) = Vo ^1 — ^ (Brane Inflation (BI)) 1671 . 

and V(tp) = Vo ^ (Dynamical Supersymmetric Inflation (DSI)) fl 68 l . 

since one has £ = and the minimum of a s is obtained at n s = 1 — ^ with 

r = 0.1384, once again, one gets 

3r 2 

a, >-—> -0.0018, (57) 

which is also incompatible with the runnig provided by Planck2013 at la C.L., 
because 2 ) > 1.6a. 

Finally, dealing with Large Field Inflation (LFI) given by the monomial poten¬ 
tial V(ip) = Vtyp p ll69l . and Radiation Gauge Inflation (RGI) 11701 . whose potential 
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(58) 


2 

is given by V(ip) = V 0 -^ 2 , one has 

a s > ~^( n s - l) 2 > -0.0018, 

where we have evaluated this quantity at n s = 0.9481, giving 2) > 1.6<r. 

Numerically we have tested all the 49 potentials that appear in table 1 of Ell : 
Using Planck2013+WP+BAO:ACDM+r+o. s data, 45 of the potentials in the list 
of ffB obtaining 2) > 1.6(7. The other four potentials, namely Loop Inflation 
f72l . Radiatively Corrected Higgs Inflation f73l . 6 exponential inflation lf74l and 
Generalized Mixed Inflation lf75l . 2) > 1.4(7. More restictive are the data provided 
Planck20 1 3+WP+high-U ACDM+r+a s (see table 5 of ifTOll ) where all of the 49 
models I) is greater than 1.6cr. 

From the results provided by Planck2013 and these analyses, one could find 
that, the realization of a sizable negative running behavior implies that inflation¬ 
ary cosmology may require multiple fields or a breakdown of slow roll condition. 
For instance, following the second path, in ll76ll a short violation of the slow roll 
condition has been considered, due to an appearance of a step-like potential from 
the perspective of string theory. Fortunately, for single field slow roll inflation, the 
new Planck2015 observational data lUTlfTZl , reduce the modulus of the running one 
order what allows the viablility of the majority of single scalar field slow roll infla¬ 
tionary models. For example, using Planck2015+TT+lowP: ACDM+r+o.,, where 
n s = 0.9667 ± 0.0066 and a s = —0.0126±g;gg|f at lcr C.L. (see table 4 of lfl2l L 
and the conservative constrain r < 0.25 (see figure 6 of fT2l ). for all the potentials 
depicted above one has 2) = l.lcr. Moreover, if one introduces lensing, then 2) 
will be reduced to be lower than lcr, and thus allowing single scalar field slow roll 
inflation. 

A final remark to end this section is in order: The same kind of results are 
obtained in inflation (resp. MBS), when one deals with an equation of state P = 
—p + (jv+i )<* P E71l for inflation (resp. P = pv+I)“P MBS) that depends on 
the number of e-folds N before the end of the inflationary (resp. quasi-matter 
domination) period. In fact, for theMBS one has 

= A , 2 ^ x (n.s - 1) , (59) 


which is perfectly compatible with current observation data, taking for example 
a = 12, /3 = — and A r = 12 (see for instance (49l ). Recall that in bouncing 
cosmologies a large number of e-folds is not required, since at the bounce point all 
parts of the Universe are already in causal contact, and also the flatness problem 
is avoided, because the contribution of the spatial curvature decreases in the con¬ 
tracting phase at the same rate as it increases in the expanding one (see for instance 
EH). However, in inflation, in order to solve the horizon and flatness problems, the 
number of e-folds before the end of inflation must be greater than 50 in most of the 
current models, and hence can hardly yield a sizable negative running following 
the relation d59l). 
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6 Reheating in MBS 


Gravitational particle production in the MBS has been recently introduced in j27l . 
The idea is the same as in inflationary models with potentials without a minimum 
(the so-called non oscillatory models): to have an efficient reheating, one needs a 
non-adiabatic transition in the expanding regime between two different phases in 
order to have enough gravitational particle creation. In inflation, there is an abrupt 
transition from a quasi de Sitter regime to a radiation-dominated one. During this 
transition, light particles are created and their energy density follows p r ~ a~ 4 . On 
the other hand, after the end of the quasi de Sitter phase, the inflaton field, namely 
4>, enters a kinetic-dominated period where the energy density of the inflaton field 
follows p 0 ~ a -6 |[79ir80)l. which means that the inflation energy density decreases 
faster than that of radiation, and thus, the Universe becomes radiation dominated 
and matches with the hot Friedmann Universe. 

For the case of the MBS, the background equation is depicted by the improved 
Friedmann equation such as in holonomy corrected LQC. The non-adiabatic trans¬ 
ition could happen in the contracting phase. In fact, a transition from matter- 
domination to an ekpyrotic phase with equation of state P = ujp where u 3> 1 
could be assumed in the contracting regime. The resulting model is dubbed as 
the matter-ekpyrotic bounce ll3Tll . Since in the ekpyrotic phase the energy dens¬ 
ity of the field follows p v ~ a -3 ^ 1+tJ ), which in the contracting phase increases 
faster than a -6 , anisotropies become negligible (note that the energy density of the 
anisotropies grows in the contracting phase as a -6 which is faster than those of ra¬ 
diation and regular matter, and thus, without an ekpyrotic transition the isotropy of 
the bounce is destroyed, which is known as the BKL instability |[29l ). Moreover, 
the energy density of the field also increases faster than that of radiation, which 
means that the field dominates the evolution of the Universe in the contracting 
phase, but when the Universe bounces, radiation eventually dominates, because in 
the expanding phase a _3 ( 1 +UJ> decreases faster than a -4 . 

To be more specific, we will study reheating via massless %-particles nearly 
conformally coupled with gravity, using the method developed in 11811 . It is well 
known that the energy density of the produced particles is related via the /3-Bogoliubov 
coefficient as follows f79l 



(60) 


where 



(61) 


and where £ = g is the nonminimal coupling constant. 

To perform the calculation we consider the simplest model of an abrupt trans- 
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ition from matter to ekpyrotic phase |jTQ, 


a(t) = 


fot t<t E 
(fpc(l + w) 2 t 2 + l) 3 ( 1 +™) for t > t E , 


(62) 


where w 1, to = t E — jjjp- t E is the time at which the transition occurs, 
and H e = H(t E ). Note that to has been chosen so that a(t) has continuous first 
derivative at the transition time t E . 

First of all, to remove ultra-violet (UV) divergences one has to assume that the 
second derivatives of the scale factor are continuous at all times. In this case, after 
integration by parts the /3-Bogoliubov coefficient becomes 


Pk = ~ 


(1-6Q 

4 k 2 


D —2ikr] 


*"(v) 

a(v) 


dr] 


(63) 


Now, for the sake of simplicity, we will assume as in [79] that the third derivative 
of the scale factor is discontinuous at the transition time r] E . Then, one has 


& * ■ 

As a consequence, the total energy density of the produced particles is 

Px = ^ ( 1 - 6 OV^, 

16 Ppi 


(64) 


(65) 


where p E = 3 H 2 E is the energy density at the transition time. 

In order to calculate the reheating temperature, namely Tr, first of all one has 
to define the reheating time t E as the time when the radiated energy density equals 
the background energy density. Since the background energy in the ekpyrotic phase 
for large values of t is given by p(t) = 3 .J> t2 , the reheating time is of the order 


to ~ , [?** _ 1 _ 

V Pe m 5 / 2 |l — 6^| H E 

and thus, the reheating temperature is of the order |28| 

Tr ~ p^ 4 (t,R) ~ y/|l — 6£\w 3 / 4 \M p i . 


( 66 ) 


(67) 


Note that we have written p E in terms of the Planck mass M p i = p l J 4 as follows 
p E = A 2 Mp t . Finally, choosing 1 — (if; |'«; 3 / 4 A ~ 10 -7 , this theoretical value 
can become consistent with the present observational bound lf27l . 
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7 Challenges and problems of the MBS 


In the previous section, we have shown that the analysis in the plane (n s , a s ) of 
the quasi-matter domination provides a range of values of the background field ip, 
whose corresponding theoretical values for the spectral index and its running tit 
well the Planck2013 data. What is important to stress is that in the quasi-matter 
domination regime the value of the background field determines its derivative Q, 
which means that the analysis in the plane (n s , a s ) provides a short curve in the 
phase space (p, p) of initial conditions, whose corresponding orbits (solutions of 
the conservation equation) depict a Universe compatible with the data of the spec¬ 
tral index and its running provided by the Planck experiment. Consequently, for all 
these orbits one has to calculate (numerically) the amount of non-gaussianity and 
to compare it with the Planck data lfl5l IT6 1 in order to check which of these orbits 
are viable. This is a very complicated task as we will explain here. 

7.1 Non-gaussianities in the MBS 

In ifTTIl the authors study non-gaussianities in the context of the MBS for a single 
scalar field using the formalism developed in ll82l . In that study, the background 
solution depicts a matter dominated Universe in the whole contracting phase in the 
frame of GR. As usual, a nonsingular bounce may be achieved by including mod¬ 
ified gravity effects or introducing nonconventional operators of the scalar field at 
high energy scales and hence it is expected that the nonlinear perturbations would 
only become important in the UV regime of bouncing cosmologies, which could be 
outside of the present observational window. For instance, a nonlinear treatment of 
primordial curvature perturbations passing through a nonsingular bouncing phase 
were numerically computed based on a particular model involving auxiliary ghost 
fields in ll83Tl . of which the result reveal that there is no evidence pointing to a 
manifest generation of nonlinearities in the infrared regime. The difficulties of 
calculating the 3-point correlation function during the bouncing phase were re¬ 
cently pointed out in lfl 8 l by studying a specific bouncing model: a background 
FLRW geometry with positive spatial curvature in the context of GR. One ought 
to be aware that this model differs from usual matter bouncing cosmologies since 
a quasi-matter domination in the contracting phase with positive spatial curvature 
does not give rise to a nearly scale invariant power spectrum. 

At present, it remains an open issue to quantitatively calculate the amount 
of primordial non-gaussianities produced during the bouncing phase in a generic 
bounce model. Moreover, if one goes beyond GR with the flat FLRW geometry in 
order to achieve a bouncing background, for example using the holonomy correc¬ 
ted LQC or teleparallel F(T) gravity, one has to use the very complicated second 
order perturbation equations in those theories, and try to calculate the 3-point func¬ 
tion for all the background orbits that are allowed by the (n s , ce s ) analysis. 
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7.2 Problems with particular bounce models 


Bounce models that depict the MBS with a single scalar field in LQC may also 
suffer from a gradient instability issue. For instance, dealing with holonomy cor¬ 
rected LQC (see J84l for a recent review about non-perturbative LQC), the square 
of the sound speed is given by (see for instance l48l ) 


c 


2 

s 



( 68 ) 


which means that in the bouncing phase (i.e., when the energy density satisfies 
p c /2 < p < p c ), the square of the sound speed becomes negative. In spite of 
the fact that in holonomy corrected LQC, in order to calculate the power spec- 


are 


trum, only modes that satisfy the long-wavelength condition \cjk 2 \ < 
used, it is important to realize that, in the super-inflationary regime, the UV modes 
(modes that satisfy the condition \c 2 k 2 \ — ) would suffer the Jeans instability. 

Moreover, in this super-inflationary regime, it is questionable the use of the linear 
perturbation equations, what could invalidate the results obtained about the value 
of the power spectrum because, as we have shown in formula (361, it is calculated 
throughout the whole background regime. This issue is also associated with the 
unclear trans-Planckian physics that is expected to be described by the full theory 
of LQC I l86ll85l . However, one should be aware that these dangerous modes, since 
they are in the extremely UV scales, do not appear in the observable window of 
today’s experiments. 

Dealing with F(T ) and F(R) gravity, the square of the sound speed is always 
positive. Thus, one can see that the issue of gradient instability is model depend¬ 
ent. However, as we have already explained in section 3, teleparallelism suffers 
from the problem that the main invariant, the scalar torsion T provided by the 
Weitzenbock connection, is not at all an invariant, in the sense that it depends on the 
choice of the orthonormal basis in the tangent bundle (the tetrad). More precisely, a 
local Lorentz transformation in the tangent bundle applied to the orthonormal basis 
changes the value of the scalar torsion. In contrast, the scalar curvature R provided 
by the Levi-Civita connection is a true invariant, but the problem with modified 
F(R) gravity (recall that GR, i.e. F(R) = R, forbids bounces for geometries with 
a flat spatial curvature) is that the current very complicated bouncing models, ob¬ 
tained using the reconstruction technique lf87l , do not support a matter domination 
in the contracting phase, and thus the power spectrum is not scale invariant lf88l . 

Another different possibility is to consider the MBS, for flat spatial geomet¬ 
ries, in the context of GR by including extra matter fields, which can violate the 
Null Energy condition, such as a ghost condensate field ll89ll or a Galilean type one 
lf90l l911. In these models, one ought to be aware of the potential graceful exit prob¬ 
lem as well as the gradient instability issue. For example, in a concrete cosmology 
of two-field matter bounce j92|, it was found that the model can be free from these 
dangerous issues with cosmological perturbations evolving through the nonsingu¬ 
lar bouncing phase almost unchanged li93t . Accordingly, the tensor-to-scalar ratio 
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could be too large to explain the CMB observations and may require a curvaton 
mechanism Ii23il to give rise to an enhancement on curvature perturbations from 
entropy fluctuations. 

8 Conclusions 

In this article we provided the main requirements that ought to be satisfied in a the¬ 
oretical model of the MBS so that it could be a viable alternative to the inflationary 
paradigm. These requirements are from both observational and theoretical aspects. 
In particular, we performed a detailed analysis of the primordial perturbations in 
an extended scenario of the MBS, in which the equation of state parameter slightly 
deviates from zero in a dynamical way. Due to this dynamical deviation, we find 
that the model can naturally yield a non-vanishing spectral index n s and running 
parameter a s . As a result, the parameter space of the MBS models is greatly en¬ 
riched. 

We applied this extended scenario into the theories of holonomy corrected LQC 
or teleparallel F(T) gravity and then examined the values of tensor-to-scalar ratio 
r, as well as studied the gravitational reheating process. By combining the observa¬ 
tional bounds upon the perturbation parameters A s ,n s , a s and r, we conclude that 
the present MBS models can be constrained very well. In particular, we explored 
the possibility of realizing a sizable negative running in bouncing cosmologies, 
which is expected to be tested in latest observational data (namely, the Planck data 
released in 2015 (Til). Additionally, we also commented on several unresolved 
issues including nonlinear perturbations during the bouncing phase, the Jeans in¬ 
stability arising in the UV regime, and the Lorentz invariant in F{T). These issues 
are not only important in constructing specific models of bouncing cosmologies, 
but they also put forward theoretical challenges to understanding quantum gravity 
theories. 
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